Entanglement of formation from optimal decomposition 
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We present a new method of analytically deriving the entanglement of formation of the bipartite 
mixed state. The method realizes the optimal decomposition families of states. Our method can lead 
to many new results concerning entanglement of formation, its additivity and entanglement cost. 
We illustrate it by investigating the two-qubit state, the separable state, the maximally correlated 
state, the isotropic state and the Werner state. 
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The rapid progress in quantum information theory 
(QIT) has made the entanglement the key resource in 
quantum mechanics. One of the most fundamental ques- 
tions in QIT is how to quantitatively describe the entan- 
glement by entanglement measures 

The main proposed entanglement measures so far in- 
clude the entanglement of formation (EOF) Ef (and the 
corresponding asymptotic generalization, the entangle- 
ment cost Ec) and distillable entanglement Ed 2]. They 
respectively present the entanglement of a state in terms 
of EPR pairs that is required to create it and that can 
be extracted from it. All the measures are equivalent 
to the von Neumann entropy E{\ip)) = S'(Tr^|V')('0|), 
S{p) = — Trp logp for pure bipartite states [3]. 

The analytical calculation of entanglement measures 
is important. The entanglement measure quantifies the 
quantum correlation in the state used for any quantum- 
information task, e.g., an EPR pair is both sufficient 
and necessary for quantum teleporation. A straight way 
to prove an important conjecture that the formation of 
mixed entanglement is irreversible 0, 0] j is to calculate 
the entanglement cost and distillable entanglement. If 
the conjecture is true, it will mean that the mixed en- 
tanglement manipulation is irreversible under local op- 
erations and classical communications (LOCC) Q. On 
the experimental side, there has been recent methods of 
estimating entanglement measures like the EOF [7] . En- 
tanglement measures have been also shown to be useful 
in many other fields such as quantum phase transition, 
spin-boson model and optical lattices 3 . 

Despite the importance of entanglement measures, it 
is very difhcult to derive their analytical expressions for 
mixed bipartite states by reason of the involved optimiza- 
tion problem. Main progress was made in the derivation 
of EOF 0, [nL There are also a few results on en- 
tanglement cost 12, 14 1 and distillable entanglement 



In this Letter we derive the EOF by means of the op- 
timal decomposition (OD) of states. We show that a 
known OD from a state always leads to an OD family 
containing infinite number of states, and the EOF of any 
state in it is computable. Next, new OD families can be 
derived from the tensor product of OD families with ad- 
ditive EOF, so we can continue to generate states with 
computable EOF. Differing with the former skills, our 



method mostly computes the EOF of states which are 
not symmetric and whose subsystems have different di- 
mensions. It greatly enlarges the family of states whose 
EOF is computable and additive. Moreover, Our method 
can evaluate the entanglement cost for many states, while 
the latter is important but very difhcult to compute so 
far. The new results on additivity of EOF help get insight 
into the classical capacity of quantum channels [l^. 

We give the method of finding out the OD of a state 
with known EOF. We determine the ODs of rank-2 two- 
qubit states with additive EOF and two important states 
in QIT, the Werner state [l7| and the isotropic state 
0. We describe how the ODs of these states lead to 
computable EOF and entanglement cost when the EOF 
is additive. We propose the schemes for experimentally 
generating the Werner state and checking the EOF. The 
separable states and maximally correlated (MC) states 
of arbitrary dimension, whose EOF turns out to be ad- 
ditive, are also analyzed. They provide strong evidence 
that the entanglement cost is strictly larger than the dis- 
tillable entanglement. 

Let us start by recalling the definition of EOF. A mixed 
state p — X]iPilV'i)('0i| can be prepared in many ways of 
ensembles of pure states {pi, The EOF for p is 

defined in the way such that 
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Suppose the ensemble {qi, \(j)i),\/qi > 0} realize the EOF 
and it thus constitutes the OD of p. By randomly chang- 
ing the probabilities, we refer to the set of states £p 
constituted by the ensemble {r^, Vr^ > 0} as the 
OD family of p. That is, a state ct S £p if and only 
if (7 = '^iri\(j)i){(j)i\,ri > 0. For simplicity, we denote 
£p = or ep = {ro,ri, |(/)o), |0i), ■■•} through- 

out the paper. The EOF of any state in an OD family is 
computable [llj, namely Ef{a) = X^i ''i-^d'/'i))- addi- 
tion, the subset of an OD family is also an OD family. 

On the other hand, the OD family is closely connected 
to the additivity of EOF of p. The latter means for any 
bipartite state a, it holds that Ef{p®a) = Ef{p)+Ef{a). 
The problem that whether EOF is additive is important 
[l^ . for it will lead to the equation Ef = Ec- The follow- 
ing facts easily follow from the definition of additivity. 

Lemma 1. (1) The EOF of any state in the OD family 
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of p is additive if the EOF of p is additive; (2) If the EOF 
of two states are additive, then the EOF of their tensor 
product is additive. 

Proof. (1) Let the OD family of phe Sp — {r^, 
We are always able to write p = '^^PiPi, with each state 
Pi G Ep and \fpi > 0. For any state a we get 

Ef{p) + Ef{a) = J2P^{Ef(p.,) + Ef{a)) 

i 

> ^piE}{p,®cj)>E}{p®(T),{2) 

i 

where we have used the fact that the EOF is subadditive 
and convex [l| in the inequalities. Because Ef{p) is ad- 
ditive and the state pi is arbitrary, we conclude that the 
EOF of any state in Ep is also additive. 
(2) Suppose Ef{pi) and Ef{p2) are additive. For any 
state cr, it holds that Ej{pi®p2®<7) = Ef{pi)-\-Ef{p2® 
a) = Ef{pi® p2) + Ef{a). So the EOF of the tensor state 
pi ® p2 is additive. This completes the proof. 

Proposition 1. Let Ep^ = {pi, and Ep^ = {qi, 
be two OD families and at least one of Ef{pi) and Ef{p2) 
be additive. Then the set e^ = {^ij A'^i'l'j)} is an OD 
family. If both of Ef{pi) and Ef{p2) are additive, then 
the EOF of any state a' S Ea is additive. 

Proof. Suppose the ODs of the states pi and p2 are 
constituted by the ensembles {p^, \ipi)} and {g-, \4>i)}, re- 
spectively. The OD of the tensor state pi (g) p2 is then 
constituted by the ensembles {p'iQj , \'ipi4>j)} because of 
Efifi (g) P2) — Ef{pi) + Ef{p2). This makes the set 
Ea = £pi^p2 = {'TijM'i't'j)} an OD family, too. When 
Ef{pi) and Ef{p2) are additive, equivalently Ef{pi(^ P2) 
is additive by (2) in lemma 1, the fact that the EOF of 
any state ct' G e^ is additive is immediately derived from 
(1) in lemma 1. I 

Proposition 1 is very useful for entanglement measures. 
In the following text the known EOF of the states p^'s 
on the spaces Ha^ (X" Hb^ 's respectively, helps generate 
the OD families Ep^ = {py , |0ij)}'s such that the EOF 
of any state ai G Ep. is computable; that is, Ef{ai) = 
Pij-^i\4'ij))- Our method indeed generates a family of 
states with computable EOF by using only a state with 
known EOF. When Ef{ai) is additive for some state ai, 
by proposition 1 we can generate the new OD families 
Ecr,<s,cr, = {Pi]k, \4>ik)\4>jk)} on the space HA^Aj Hb^Bj 
and the EOF of any state ct^ G e<Ti(g)<Tj similarly reads 
Efi^ij) = J2kPi3kEi\4'ik)\4'jk))- In the same vein, we 
can generate more OD families Ea-i^a-j^-.-^a-k on the space 
HAiAj...Ak HBiBj...Bk and calculate the EOF of states 
in it. When Ef{ai),Ef{aj), ...,Ef{ak) are additive, by 
proposition 1 they will lead to new states whose EOF is 
additive and thus the computable entanglement cost. 

For example, we consider the two-qubit MC state 

Pp,e =p\ipe){ipe\ + (1 - p)|'07r/2-e)(V'7r/2-el, (3) 

where \^pg) = cos 6*100) + sin 6*1 11) and |0), |1), ... are com- 
putational basis. It follows from [iJI that Ef{pp^g) = 



EciPp.e) = h{cos'^9) with the function h{x) being the 
binary entropy function. So the set Ep^ g = {p, 1 — 
p, li'e), \ipTT/2-e)} is an OD family. Thus by proposition 
1 we can compute the EOF and entanglement cost for 
any state p in the OD family £pp^,e^0Pp2.f>2'»-'»Ppr..o„ = 
{pi,P2, ...,P2",\'4^e,'4'e2--ipeJ, \il^^/2-9i'4^62-i'eJ, ■-, 
lV'7r/2-eiV'77/2-e2---V'7r/2-e„)}, namely Ef{p) = Ec{p) = 
S"=i h^cos^Oi). Remarkably, the EOF of different states 
in the OD family do not change with the probabil- 
ity distributions. So we have proposed a new family 
of MC states and thus Ef{p) is additive. As the dis- 
tillable entanglement of the MC state is computable 
[3], we thus obtain that the difference between the 
entanglement cost and distillable entanglement of p G 
«^Ppi.ei«'PP2,£)2®-«'Pp„,e„ on the space Ha «> Hb is 

n 

E,(p) - Ea{p) = J2 Kcos^S,) - S{Ttap) + S{p). (4) 

i=l 

This expression indeed represents the entanglement that 
cannot be distilled from the state p under LOCC 0]- 
When the conjecture that Eqs. (4) is strictly larger than 
zero is true, it will mean the asymptotic entanglement 
manipulation of p is irreversible, which is essentially dif- 
ferent from the case of pure states Much numerical 
calculation has been done and they supported the con- 
jecture. Furthermore, it turns out to be true by plotting 
Eqs. (4), when the number of parameters in p is less. 

Similarly one can verify a more general two- 
qubit OD family e^ — {q,l ~ q, y/pcosd\00) + 
Vl~^sin6l(x|00) + y|01) + z|ll)), y^sin 6l|00) - 
VI -pcos6'(x|00) -I- ylOl) + z|ll))} up to the normal- 
ization factors, p,x,y,z G [0,1], -\- -\- = 1 and 

tan 6 = — — ^ . The EOF of any 

state in Ea is equal to h{\ -I- ^■\/l — 4(1 — p^x^z^') 
and it is additive [l^ . As the state pp,d in 
Eqs. (3) belongs to e^, the newly constituted 
OD family Ea^®a2®---®<j.„, provides more states p 
whose EOF is computable and additive, namely 
Es{p) = E,{p) = Y:UK\ + iVl-4(l-p.)X'^f). 
Generally, 9] has given a programmable way to find out 
the OD of any two-qubit state p. It provides abundant 
resource for the generation of new OD families, in which 
the EOF of the states is computable by proposition 1 . 

In this case, we can obtain more states with com- 
putable EOF by mixing Ep with other OD families, e.g., 
the family Es consisting of all separable states ps. As the 
entanglement cannot be increased under LOCC and thus 
Ef{p^ Ps) > Ef{p), it then follows from the subadditiv- 
ity of EOF that 

Lemma 2. The EOF of any state p G Eg is additive. 
The set of separable states Es is an OD family. ■ 

Although the separable state is classically established 
[l3| , we emphasize that lemma 2 is by no means a trivial 
result. We consider the 4x4 state a in the OD fam- 
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ily {p,l-p,|00)|V'fl),MK)mi>|^^/2_,)} 0, ^hich 

is a subset of the OD family Ep^^p^g- So we get that 
Ef{a) = Ec{<7) = h{cos^9) and it is additive. However, 
these results cannot be derived from any existing theory 
since the partial trace over a system of a does not equal 
to an entanglement-breaking channel [isj . In addition, 
the tag states |00) and l^^^l^^^ cannot be explic- 
itly distinguished iJLSi- Generally given an OD family 
Ep — {pi, by proposition 1 one can always consti- 

tute the new OD family Ep^^p — {pi, \(j>i)\ipi)} with each 
state \(f)i) in randomly product form. The EOF of any 
state a € Sp^^p thus reads Ef{a) — J2iPi^{\^i))- This 
value will coincide with Ec{(j) when Ef{p) is additive. 

We have known that the OD of the two-qubit state 
p is constituted by pure states with identical amount 
of entanglement [9|. However it is not the case 
for states of high dimensions, as the following ex- 
ample illustrates. Suppose = X^itTo^ '^d**) J'/') — 
(Ef="o cf )-'/' . COS e = (ZUjy^i^ = 

(sr-d-1 2\l/2 ,,_„,_ /l I -l+p+pcos2e _ 

^10 = \l-2- ;TT°l% ^'^^^ > 0, d > /, we have 

Lemma 3. Consider the MC state p = p\4'){i'\ + 
(1 - p)\(t)){(i)\. The OD and EOF of p respec- 
tivel y read Y.]=Q{uM^\i)) + un^l - p\(f))){uiQ^{ij:\ + 
Mil VI - and 

d-l /-I 2 2 

Ef{p) = -pE c| logcf - (1 - p) E ^ log ^ 

1=0 i=0 

+Hl + 5^1 -P^sin" 29)~ph{cos^ 9). (5) 
Proof. It follows from the definition of EOF 9] that 

Ef{p) ^ mm 2_^q^E[ ^ — , (6) 

where [wy] — [t/o, wi] are the first two columns of a unitary 
matrix and ^Jql is the normalization factor. Notice every 
decomposition of p is expressed in Eqs. (6) owing to the 
linear independence of I '0) and 10). Reduction of Eqs. (6) 
leads to the optimization problem equivalent to that for 
the two-qubit MC state, whose entanglement has been 
derived by @. One can check that the entanglement of 
the proposed OD of p coincides with Ef{p). ■ 
Lemma 3 actually provides a new family of MC states 
with computable EOF. One can easily check that the 
amounts of entanglement contained in the two pure states 
constituting Sp are always not equal, which totally differs 
with the OD families of two-qubit states and separable 
states. It thus by proposition 1 generates new OD family 
£pi0P20---0p„ in which the pure states have different en- 
tanglement, correspondingly the states in Epi0p2<s---0p„ 
are not equally entangled when the probability distribu- 
tions change. So this OD family generates more complex 



states with computable entanglement cost. Of course 
more OD families can be produced by using the results 
on two-qubit states and separable states. 

Next, we show that there is indeed a gap between the 
distillable entanglement and entanglement cost for p. It 
follows from [l^ and Eqs. (5) that 

E,{p) - Edip) = + iVl-p' sin^ 29) ~ h{psm^ 9) 
+h{\ + \^/l- 4psin^ 9 + V sin^ 9) , (7) 

which turns out to be strictly larger than zero unless p = 
or 1. A simple approach to this inequality is by plotting 
Eqs. (7). Differing with the qualitative derivation (e.g., 
[J]), our result analytically shows that one can distill the 
same entanglement from p as Ec{p) if and only if p is 
pure. The undistillable entanglement is explicitly given 
by Eqs. (7). It strongly supports the conjecture that 
entanglement distillation is irreversible in general 

Another proof for OD family containing different en- 
tanglement concerns an important state in QIT, the 
isotropic state pp = + -^3rlV'+)(V'+M'/'+) 

Y^f^i [IBl- It turned out that for pp with d — 3 
and F > 8/9, Ep^ consists of a maximally entangled state 
and nine states p obtained by twirling and they are thus 
equivalently entangled, i.e., E{p) — —1/3 + log 3 fld^ . 
We prove that Sp^ of high dimension d,F > {Ad — 4)/(i^ 
is similarly composed of \ip^) and many states p ob- 
tained by twirling, each of which contains entanglement 
E{p) ~ ^^\og{d — 1) -I- log d. This conclusion follows 
from the OD oi pp such that 

l^.>^^Ea.N>Ear.N> + ^ER. (8) 

^ l—l l—l ^ l—l 

We demonstrate the L x d coefficient matrix [an] in 
the case of odd d > 3. We regard the n x 1 vectors 
dj,j — 1, ...,d as the nonzero entries in [an], each "row" 
of which consists of d + 1/2 nonzero entries and d — 1/2 
zero. The subscript of dj marks its column in [an]. 
The first row of [an] consists of [di, d^+x/2i 0, 0], 
the second [di, dd_i/2, 0, 0^+3/2, 0, 0], and the 
last [0, ...0, d(i+i/2, dd]. Counting all kinds of combi- 
nations, there are in all {j^_^!^i^2} ^^ws and hence L = 

"(d+'i/2)- Each dj has the form (y^^, \f^ef'^\ 

e-''3^("~i)*)^, where the natural numbers //s are 

required to be + fj — fk — fi if and only if i = fc, j = Z 
or i = l,j = k. It can be done, e.g., by choosing 
fi — m^,m > l,n > 2m'^ — 4. One can deal with the 
case of even n similarly and verify Eqs. (8). So we have 
given the method of generating numerous Epp , which is a 
remarkable character of the isotropic state. It makes Ep^ 
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a much stronger OD family compared to the former re- 
suhs, when it is used to create the states with computable 
EOF by proposition 1. Specially, the EOF of any state 
in Epp has the form p(2^1og(d— 1) + log d) + (1 — p)log d 
with the probability p £ [0,1]. In addition, the OD of 
Pf implies that one can generate the isotropic states 
by classically mixing the ensemble of MC states p — 
p\ipi){tpi \ + (1 — p)\tjj~^){tp'^\ obtained by twirling. 

On the other hand there exist the states of high dimen- 
sions whose OD consists of equivalently entangled states, 
e.g., the Werner state p^,, which is important in QIT and 
has been extensively investigated Ull . A dxd Werner 

state has the form p^ = + ^jT^J^uto Ni)0'«l- 

The EOF of an entangled Werner state {F e [-1, 0)) has 
the analytical expression Ef[pw) = h{l/2 + \/T~~F^/2) 
From it we can derive the OD of p^j such that 



d-l 3 

P^=2d^ T, T, \^i3k){'4'^jk\, 
i>j=Q k=0 



+2"fciyg||i7>+2ufc2/ 



1-F 
2 72 



(rf-l)(f +1) + 

2d+2 72 ' 



2d+2 U-J 

/ -1/2 1/2 1/2 1/2 \ 
1/2 -1/2 1/2 1/2 
i/2 i/2 -1/2 1/2 
y i/2 i/2 1/2 -1/2 / 



(9) 



One can easily verify the proposed OD which indeed gives 
rise to the Werner- OD family — {pijk, IV'iife)}; s-nd 
each pure state \ipijk) has identical amount of entangle- 
ment Ef{pw)- It thus implies that the EOF of any state 
in Ew is equal to Ef[pw), no matter how the probabil- 
ity distributions change. This is similar to the case of 
two-qubit states. We can use the Werner-OD family to 
construct more OD families with computable EOF by 
proposition 1. It is a new function of the Werner state. 
Moreover, one can check that any state J2ijPijPij ^ 
is negative partial transpose (NPT). This helps infer the 
EOF and irreversibility of NPT bound entangled states, 
if it really exists 22l |. 

We notice that for each pair of fixed i,j, the state 
Pij € £w is just a 2 X 2 Werner state p^o- Such an 
OD is interesting in the sense that we can experimen- 
tally prepare a, d x d Werner state by means of clas- 
sically mixing many states Pwo^s up to unitary opera- 
tions with identical probabilities pij 's, equivalently mak- 
ing a state p^a go through the unital channel A{p) = 

d^J2i>j=o'"i3 '^^jPwovlj 'S) wlj. Here each pair of 
unitary operations Vij and Wij acts on a 2 x 2 space, so 
they can be indeed regarded as Pauli operations / and ax ■ 
As the state p^o has been reahzable [21[, the proposed 
experiment is probably realizable by current techniques. 
More importantly, the experiment will verify that the 
Ef{pw) is indeed the minimal entanglement required to 
create a Werner state of high dimension. 



To summarize, we have presented the OD method of 
deriving the EOF, additivity and entanglement cost for 
many states. Our method is flexible and could yield more 
results on entanglement measures. It also helps generate 
the Werner state and check the EOF experimentally. 
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